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Abstract 

Following the paper by M. Combescure [Ann. Phys. (NY) 204, 113 (1990)], we apply the 
quantum singular time dependent oscillator model to describe the relative one dimensional 
motion of two ions in a trap. We argue that the model can be justified for low energy excited 
states with the quantum numbers n <C Umax ~ 100, provided that the dimensionless con- 
stant characterizing the strength of the repulsive potential is large enough, ~ 10^. Time 
dependent Gaussian-like wave packets generalizing odd coherent states of the harmonic oscil- 
lator, and excitation number eigenstates are constructed. We show that the relative motion 
of the ions, in contradistinction to its center of mass counterpart, is extremely sensitive to 
the time dependence of the binding harmonic potential, since the large value of results in 
a significant amplification of the transition probabilities between energy eigenstate even for 
slow time variations of the frequency. 
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1 Introduction 



A quantum system described by the singular oscillator Hamiltonian 



is one of rare examples admitting exact solutions of the Schrodinger equation. In the case of 
a constant frequency co, the eigenfunctions can be found in textbooks [0. Generalizations of 
these eigenfunctions to the case of a time dependent frequency were obtained for the first time 
in IQ. Using these solutions, an exact propagator of the Schrodinger equation with Hamiltonian 
(H) was found in An elegant method of calculating propagators of time dependent quantum 
systems, including (||), on the basis of quantum operator integrals of motion, was developed in 
[^. The same propagator was obtained in the framework of the path integrals approach in |^] (a 
constant frequency) and (a time dependent frequency). A detailed analysis of solutions found 
in Q was performed in Ref. 0, where two families of time dependent wave packets generalizing 
coherent states of the harmonic oscillator were also constructed. An asymmetrical singular 
nonstationary oscillator in two and three dimensions was considered in P] (see also |10| for the 



stationary case). A systematic study of wave packets and coherent states in the centrifugal and 
other anharmonic potentials was given in |11]. 



The exact solvability of the Schrodinger equation with Hamiltonian (||) is explained by the 
fact that this Hamiltonian is a linear combination of three generators of the su{l, 1) algebra, 



2 

— + ^, X , xp + px. (2) 
2fi x^ 

The propagator and solutions of the Schrodinger equation with a generic Hamiltonian, which 
is a linear combination of generators (|2|) with arbitrary time dependent coefficients, were found 
in (see also U], where other representations of the sm(1,1) algebra were considered). 
The universal invariants and generalized uncertainty relations for this generic case were studied 
in [|l5| . A special case of periodic time dependent coefficients, including the adiabatic (Berry's 
phase) and delta-kicked limits, was considered in Recently, many results obtained in the 
cited papers were rediscovered, see e.g. |l8|, p^ . 

A distinguished role of Hamiltonian (jl|) is explained by the fact that, in a sense, it belongs 
to a boundary between linear and nonlinear problems of classical and quantum mechanics. For 
this reason, it was used in many applications. For example, it served as an initial point in 



constructing interesting exactly solvable models of interacting A^-body systems [20[. It was used 



also for modelling diatomic [12| and polyatomic |21] molecules, and for the evaluation of the 
accuracy of the quasiclassical approximation for nonquadratic Hamiltonians [^]. Studying the 
limiting case g = of the generalized coherent states of a singular oscillator led to introducing a 
concept of even and odd coherent states which are considered nowadays as simplest models 
of quantum macroscopic superpositions — the so called Schrodinger cat states [p3||. 



Even and odd coherent states of photons were created in cavity experiments [^]. Another 
superposition of coherent states of the Schrodinger cat type, similar to even and odd coherent 



states, can be created from an initial coherent state in a nonlinear Kerr medium |25]. This su 



perposition has the same Poisson quantum number distribution function as the usual coherent 



2 



state, so it is "more classical" |26] than even and odd coherent states, which have highly oscillat- 
ing quantum number distribution functions (strong oscillations of photon distribution functions, 
like in squeezed states are believed to be a characteristics of nonclassical states |]29|| ). 

The photon statistics of parametrically excited even and odd coherent states (which are super- 
positions of squeezed states) was discussed in |3C]. A quantum sling mechanism of constructing 
squeezed and correlated states from usual coherent states was suggested in |31]. 

Unfortunately, photon nonclassical states are very sensitive to the interaction with an en- 
vironment, and decoherence phenomena may destroy the specific interference pattern of the 
Schrodinger cat superposition states very rapidly (for example, an increase of temperature leads 
to a decrease of oscillations of the photon distribution function in squeezed and correlated states 
[^2| ). Much more stable even and odd coherent states can be created for trapped ions (for exam- 
ple, in the Paul trap P3]), as was demonstrated both theoretically p3, p5[ and experimentally 



[ p6| . Different aspects of quantum motion of ions in traps were discussed, e.g. in |26, In 
particular, various nonlinear modifications of coherent states []39| could be created for trapped 
ions ||4^. 

Here we wish to show that the limiting case of very large values of the dimensionless param- 
eter ^ig/fi^ is also interesting, because it may have some relation to the problem of the relative 
motion of ions in electromagnetic traps. Indeed, the ions in the Paul or Penning traps move 
in time dependent harmonic potentials, and the stable configurations correspond under certain 
conditions to an effectively unidimensional motion along the trap axis |41]. In such a case, a 
relative motion of two ions can be described as the motion of a single particle (with the reduced 
mass) in the Coulomb repulsive potential V{x) = e^/x. Combescure [H] proposed to replace 
the real Coulomb potential by the singular oscillator potential Vef = g/x^ . Using the solutions 
of the Schrodinger equation in the same form as in |Q, ^ and considering the harmonic time 
dependence of the frequency uj{t), he made some qualitative conclusions concerning the relative 
motion of the ions in the Paul trap. However, he did not give any estimations which could 
justify his model. The aim of our article is to provide such a justification for large values of 
parameter g, and to consider specific features of solutions in this special case, stressing quali- 
tative differences between the relative and center of mass motions. We shall demonstrate that 
the transition probabilities between energy eigenstates of the relative motion are significantly 
amplified in comparison with that of the center of mass motion, even when time variations of the 
frequency are very small. This result may be interpreted as an indication to a high sensitivity 
of the relative ion motion to small time variations of the binding trap potential, which could 
result in some kind of a chaotic behaviour in the case of the real Coulomb interaction potential 



U 43, Note that an irregular (although not chaotic, in the strict sense of this notion) be- 



haviour of time dependent linear quantum systems described by su{l, 1) invariant Hamiltonians 
with (quasi)periodic 5-kicked coefficients was studied in [45, E^] (see also 



2 Justification of the model 



Let us compare two potentials (// is the reduced mass of two ions). 



V{x) = }-fiu^x^ + -, (3) 
2 X 
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V,{x) = \i.u?^x' + A. (4) 

If the energy of a particle is not too large (as in real traps, where ions are cooled), then it 
performs small oscillations near the equilibrium coordinates satisfying the equations V' [x) = 
and Vg{x) = 0, namely, 

The frequencies of oscillations are determined by the second derivatives of the potentials in the 
equilibrium points. 



Demanding the equilibrium positions and the effective frequencies of the real potential, @, and 
its simulation, (Q), to coincide , we obtain the following parameters of the effective potential 
describing the relative ion motion: 



x/3 _ 3 

The minimum is twice less than that of the real potential 



1/3 

= 5 = i ( I • (5) 



but this shift seems not important. 

The solutions of the Schrodinger equation depend on the dimensionless parameter 

g* = 2fig/n\ (6) 

Taking the rc-dependent part of the trap potential energy as Vt{x) = eC/x^/(2L^), U being the 
voltage applied to the cap electrodes separated by the distance 2L, we have = eU / {2fiL?') 
(since /x = mj/2, where rrii is the ion mass). Therefore, 





"4Ry 




1/3 


j 4m 


eU 







(7) 



where m is the electron mass, Ry is the Rydberg constant, and a_B is Bohr's radius. For typical 
trap parameters, /x/m ~ 10^, U ^ 100 V, L ~ 1 mm, we have g^, ^ 10^'^. 

A specific feature of the singular oscillator potential is the equidistant character of the energy 
spectrum. Since the real interaction potential does not possess such a property, the replacement 
of (^) by (^ can be justified, provided that the corrections to the energy spectrum due to 
the anharmonicity of the potential (^) are small, comparing with the effective oscillator energy 
difference hQf, = \/^hLO. The cubic anharmonic correction to the potential (^) in the vicinity 
of the equilibrium point is — (e^/Xg) {6x)'^, 5x being a small deviation from the equilibrium 
position Xe- Assuming ;ur2g(5x)^ = nhQe for the nth excited effective harmonic oscillator energy 
level, one can verify that the anharmonic corrections can be neglected under the condition 



n <C n^ax ~ 3 ( ^3- j '^Sgl^'^. (8) 



For the same trap parameters as above, Umax ~ 100. Consequently, the singular oscillator model 
can be applied for describing the low energy excited states of the relative ion motion. 
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3 Time dependent solutions for a large repulsive constant 



A complete set of orthonormalized solutions to the Schrodinger equation with Hamiltonian (||) 
can be chosen in the form [Q, ^ 



d+l 



n! 



he^J T{d + n + l) 



— exp 



fix 



(9) 



where d = ^^/T+'^g^, and L'^{z) is an associated Laguerre polynomial. (Hereafter we write 
a; instead of ujg, assuming that the frequency rescaling has been done.) The time dependent 
complex function e(t) is a special solution to the classical harmonic oscillator equation, 



e + Lo'^{t)e = 0, 



satisfying the normalization condition 



e e 



2i. 



(10) 



(11) 



If w = const, the solutions (|9|) go to the energy eigenstates, provided that function e{t) is chosen 
as e(t) = o;^^/^ exp (iut). 

Expression (^) is written for x > 0. Its continuation to the negative semiaxis (since x 
is the relative coordinate of two ions along the trap symmetry axis, it can be both positive 
and negative) depends on the symmetry properties of the two-particle wavefunction determined 
by the kind of statistics (Bose or Fermi, if the ions are identical), so that one should write 
^(x) = ib\I'(— x) for X < 0. There are no problems with the continuity of the solutions or 
their derivatives at x = 0, provided that d > 1/2 (i.e. g > 0). Extending the wave function 
to the whole axis, one must only change the normalization by the evident factor ^/2. Authors 
of the recent paper |48] claimed that there is a nontrivial influence of the parameter d on the 
possibility of existence of the bound states of two ions in the Paul trap, depending on the type 
of quantum statistics. However, this conclusion was a consequence of incorrect manipulations 
with the expressions for the wave function at x > and x < 0, which resulted in the coordinate 
independent factor libcos((i+l/2)7r in the wave function. This factor turns into zero for specific 
values of parameter d, which was interpreted by the authors of |Q as a disappearance of bound 
states. They forgot, however, to normalize the wave function. If they did this, then the incorrect 
factor would disappear, not the bound states. 

Solution (|9|) is an eigenfunction of the time dependent integral of motion (a generalization 
of the harmonic oscillator number operator) 



B 




\e\' \ 1-^1+ ^|ep2;^ - ^Re (ee*) {xp + px) 



with an eigenvalue 2n + d + 1 . There exists another integral of motion Q , 

1 



A 



[ep — e/ixj + 2e 



X' 



(12) 



(13) 



which is a generalization of the square of the harmonic oscillator annihilation operator a?. A 
normalized eigenstate of operator (|l3|), A\a) = a^|a) (a being an arbitrary complex number). 
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reads Q 



{x\a) 



E 



a 



2n 



\2he 



nl/2 



^62l,(|a|2) 



exp 




(14) 
(15) 



where Jd{x) and /^(x) are the Bessel function and the modified Bessel function, respectively. 
Similar (time independent) eigenstates of the lowering generator of an abstract sm(1, 1) algebra 
in another specific representation were considered by Barut and Girardello |^0|. Various gener- 
alizations and concrete realizations of these states were discussed, e.g. in |5^, 54 1, where 
many references to other publications can be found. 

Taking d = 1/2 (i.e. g = 0) and replacing the reduced mass, fi, by the total mass of two 
ions, M, we arrive at the expressions for the wave functions describing the odd states of the 
center of mass motion in the harmonic potential with frequency uj (remember, however, that the 
frequency of small oscillations of the relative distance between ions is \/3 times higher than the 
frequency of the center of mass motion) . The even wave functions of the center of mass motion 
are obtained by means of a formal substitution d = —1/2. 

The time dependences of the mean values of the operators (pi) in any state are as follows. 




n 

m 

h 

2 



{xp + px) 



|e|2(^) - Re (e2(it) 
€\^{B)-Re (e2(it) 
2h Re{ee*){B) - Re (ee(it) 



In particular, the mean value of the operator B in the state ( p^ equals 

,2/^(l«P) 



{a\B\a) 



1 + a 



Idi\a\ 



(16) 
(17) 
(18) 

(19) 



Formulas ([T^) and ([1^) can be simplified for large values of coefficient d. Taking into account 
the asymptotics of the Bessel functions of a large index , 



2TTld{z) « (d^ + z^) exp 



\/(P~+^ + din 



+ Vd'^ + 



(i>l, d>z. (20) 



and calculating the power series expansion of the exponential argument in ( pO[ ) up to terms, 
we obtain 



{a\B\a) = l + d + 



\a\ 
'2d 



(21) 



Consequently, the energy difference between the states \a) and |0) becomes noticeable only for 
\a\^ ~ d. In such a case, due to Eq. ([l^ , the mean value (x^) oscillates in the vicinity of the 
average position (xq) = {h / m)\e\'^ (B) ~ 2, the amplitude of oscillations being of an order of 
"v/d ^ d. For this reason, we can simplify formula (15), replacing the Bessel function by its 
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asymptotical form following from ( pO]) and expanding the arising exponential argument up to 
the terms of the order of |a|^/(i or x^\a\^ /d. After some algebra, we obtain the expression 



(x|a) 



(27rd)- 



-1/4 



d/2 



exp 



\a\ 
'Sd 



8de2 2e^d 



where 



h 

y{x,t) = x^ |ep((i + 1). 



(22) 



(23) 



mials p^, 



In the right hand side of Eq. ( |2^ ) we recognize the generating function of the Hermite polyno- 

exp (— + 2zx^ = ^ —Hn{x). 



n=0 



Comparing it with the expansion (|r^) and using Stirling's formula, we arrive at the following 
asymptotical expression: 



^n{x,t) = (2^d) 



-1/4 ( 

V?ie22"n! 



1/2 /g*^n-fd/2 



exp 



4n2|e|4d 2^1 



The probability density in the state ( p^ is Gaussian, but with respect to variable x^: 



|(x|a)|2 = (2^d)-i/2|f^exp 



2;i2|e|4^ 



/ I 2, \ ^ 
(a x \a) = — 
m 



l + d + 



\a\ 
'2d 



Re [e*^a^) 



2dh\e\'^ J 
(24) 



(25) 
(26) 



Exact solutions to the Schrodinger equation in the form of the power-Gaussian wave packets 
were constructed in |§, @]. From the modern point of view, these packets can be considered 
as generalizations of the odd squeezed states (in the context of an abstract su{l, l)-invariant 
Hamiltonian similar states were considered, e.g. in [^5|| ; recent results and other references can 
be found in ||, ||]): 



d+l 
2\ 2 



exp (-it ) |0) 



2fix 



nT{d+i) 







z 


7^ 


d+l 


' ifi e — ze* 2 




I e — ze* 




\ exp 


2he- ze*^ . 



(27) 
(28) 



For not too large values of d, the states (15) and (^) are obviously different. However, for ^ 1 
and moderate values of x, a, and z, the function |(x|z)p is given by the same expression 
The only difference is in the formula for the mean value of x^: 



{z\x^\z) = {d+l) 



h \e — ze 
m 1 



*|2 



h 



m 



1 + 2\z\ 



Re [e^z*) 



(29) 
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For moderately excited states, parameter z must be small, since 



(z|B|z) = (d+l)i±^. 

Note that the relative width of the distributions given by Eqs. (^) and (|2^) practically does not 
depend on time, Ay/y 

max ^ 

while the position of the center of the 

distribution, Umax, niay vary essentially even with small variations of function |ep, due to the 
large value of d. As a consequence, the covering between the packets taken at different moments 
of time may be very small even for weak variations of the frequency uj{t), resulting in the effect 
of the amplification of transition probabilities, discussed in the next section. 

In Figjl] and ^ we plot the probability densities for the states ^o{x,t) and "^2{x,t) ( p^ in 
the parametric resonance case (using dimensionless units) 

9 , , 1 + k cos 2t , , , 

w2(t) = ^__, A;«l, (30) 



when an explicit analytical (approximate) expression for e{t) reads [30| 



kf kf 

eit) = cosh — e^* - i sinh — e"**. (31) 

4: 4 



We assume d = 10° and k = 0.02. 



4 Amplification of transition probabilities 

Suppose that the frequency ui{t) assumes constant values, and ujf, in the remote past and 
future, respectively (such a situation corresponds to the Penning trap, when the voltage be- 
tween electrodes exhibits some time variations), and that the singular oscillator was in the nth 
stationary state at t — > — oo. This state is given by Eq. (|9|) with ei{t) = exp^iiVit). The 

asymptotical form of the solution to Eq. (|To|) at t ^ oo reads 



— 1/2 

e{t —>■ oo) = LOr exp {iujft) — jyexp {—iojft)] , 



^ and rj being constant complex coefficients satisfying the restriction — |r/p = 1, following 
from ( [Tl| ) . The transition probability from the initial nth energy eigenstate to the final mth one 
(corresponding to frequency oOf) was found in 



^ ^Xmn\n^^n + d + l) _ \pi^rn~n\4) _ 2r)? (32) 

^mn-^ if^mn ~l~ O ~l~ 1) 



jr(i.™„ + d+l)rl™-"l(l-r)'^+i ^ , , ..„m2 



Hmn^.T {^rnn + d + 1) [|m - n| !]2 



[F {-iimn,'i^mn + d+l;\m - n\ + l]r)] . (33) 



Here fimn = min(m, n), = max(7n, n), T{z) is the Gamma function, pjf'^\z) is the Jacobi 
polynomial, F{a, b; c; z) is the Gauss hypergeometric function |^^, and r = is a reflection 



8 



coefficient from an effective potential barrier determined by the time dependence of In 
particular, the ground state excitation probability equals 



r(m + d + l) 
m\T{d+l) 



r™(l-r) 



d+l 



(34) 



In the special cases of d = ±1/2 (describing the center of mass motion in the pure harmonic 
time dependent potential), one can use the quadratic transformation of the Gauss hypergeomet- 
ric function ||49t| , 

F{a, b;a + b+ 1/2; 4z(l - z)) = F{2a, 2b;a + b+ 1/2; z), 

to reproduce the known formula for the transition probabilities between the energy levels of 
the harmonic oscillator with a time dependent frequency [^] (see also for details and the 
bibliography) , 



^ 2\t''i\{[{k + j)/2]\f 



Here k and j must have the same parity. 

Now we have to take into account that coefficient d is very large, d ~ ^Jgl ~ 10^, in the case 
of the relative ion motion, while quantum numbers m and n must be much less than 100. In 
such a case, using an approximation {b)k = b{b + 1) . . . (6 + A: — 1) 6^, 6 1, we can reduce 
the Gauss hypergeometric function 

Fia,b;c;z)^f:^-^z^ 



k=0 



{c)kk\ 



to the confluent hypergeometric function 



(ah 



(bz) 



With the same accuracy, we can write T{d + n) ~ T{d)d^. Then Eq. (|3^), together with the 
relation between the confluent hypergeometric function and the Laguerre polynomial 



<I>(-n; c + 1; z) 



{C + I)r 

lead to a simplified form of the transition probability for d ^ 1: 



' ' n 



-(r(i)l™-"l(l - rf+^ [Lj^;"l(rfi) 



(rdy 



-(1-r) 



d+l 



(36) 



(37) 



If the frequency variation is small, for instance, adiabatic, then the effective reflection coef- 
ficient r is small too (excluding the case of the parametric resonance). In this case, the Taylor 
expansion of formula (p3D yields 



w:: 



VmrX {Vmn + d + I) 



„\m—n\ 



/^mnr (//m„ + (i+ 1) [|m - n|!]2 . 



^ 2mn + {d+ \){m + n + 1) 
|m — n| + 1 



(38) 
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Consequently, if r ^ 1, and quantum numbers m, n are not too large, the transition probabilities 
of the center of mass motion (d = ±1/2) are suppressed by the factor rl"^~"L However, this may 
not be true for the relative motion, since the reflection coefficient is effectively amplified by the 
large factor d. This is seen distinctly from the special form of formula (p6D, which is valid for 
r^d < 1: 

Fig.|3| gives the dependence of on the product rd for n = 5 and m = 0, 2, 5, 10. The surface 
plots of the relative motion transition probabilities for a fixed value of d = 10^, and for 
r = 10^^10~^10-^ are given in Figs.^, |5|, and ^, respectively. Fig.^ shows that actually there 
is no transitions to other states, if rd = 0.1, and m, n < 5, while for larger values of m, n some 
structure appears in the plot. This structure becomes more complicated, when rd = 1. For 
rd = 10, the picture is quite different from the adiabatic behavior, although coefficient r is still 
very small. 

If rd ^ m, n ~ 1, then the leading term of the Laguerre polynomial in (p9|) yields 



VT- = S!:^L—{1 - rf+' ^ ^in^e-^d^ (40) 
" m\n\ m\n\ 



Eq. (|40| ) clearly shows, that for moderate values of m, n <C rd (when the model considered can 
be justified), all transition probabilities are negligibly small. This means that even small time 
variations of the binding trap harmonic potential result in the excitation of the states with high 
values of the relative motion energy, while the quantum state of the centre of mass motion is 
practically unchanged under the same conditions. This phenomenon can be interpreted as an 
indication of an instability of the relative ion motion with respect to small perturbations in the 
quantum regime. We should note, however, that the quantitative calculations of the transition 
probabilities to highly excited states of the relative motion of real ions cannot be performed in 
the framework of the singular oscillator potential model, since the condition ^ is not fulfilled 
for these states. 



5 Conclusions 

The main results of the paper can be summarized as follows. We have presented new asymptot- 
ical forms of the wave packets of the time dependent singular quantum oscillator in the case of 
a very strong repulsive centrifugal potential. We have shown that such a potential can approxi- 
mate to a certain extent the real Coulomb repulsion between two ions in a trap. We have found 
an extreme sensitivity of the transition amplitudes between energy eigenstates to small varia- 
tions of the coefficients of the Hamiltonian. It is worth noticing in this regard that the system 
of two trapped ions interacting via the true Coulomb potential exhibits under certain conditions 



a chaotic behaviour |41, 44 1. Of course, one cannot expect any chaos, in the strict sense of 
this word, in the framework of the one-dimensional exactly solvable model. Nonetheless, our 
results might be interpreted as an indication of a possible irregular behaviour (which manifests 
itself just in a high sensitivity to small changes of the initial conditions of parameters describing 
the system), if one replaces the exactly solvable potential by the realistic one (which cannot be 
treated analytically). 
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Figure 3: W^{rd) for m = 0, 2, 5, 10 ( r and d are dimensionless). 
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